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An unpaired gapless Dirac electron emergent at the surface of a strong topological insulator (STI) is 
protected by the bulk-surface correspondence and believed to be immune to backward scattering. It is less 
obvious, however, and yet to be verified explicitly whether such a gapless Dirac state is smoothly extended 
over the entire surface when the surface is composed of more than a single facet with different orientations 
in contact with one another at sharp corner edges (typically forming a steplike structure). In the realistic 
situation that we consider, the anisotropy of the sample leads to different group velocities in each of such 
facets. Here, we propose that much insight on this issue can be obtained by studying the electronic states on 
a hyperbolic surface of an STI. By explicitly constructing the surface effective Hamiltonian, we demonstrate 
that no backward scattering takes place at a concave 90° step edge. A strong renormalization of the velocity 
in the close vicinity of the step edge is also suggested. 
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A single Dirac cone emergent on the surface of a three- 
dimensional (3D) strong topological insulator (STI) is 
"topologically" protected; its existence guaranteed by 
the bulk-surface correspondence and by a nontrivial 
value of the bulk topological invariant of the strong Z2 
index. 1 ' 2 ' The existence of such a gapless Dirac cone has 
been repeatedly verified even experimentally by a num- 
ber of spin-resolved ARPES measurements performed in 
different realizations of the 3D STI 3 ~ 5 ' and is now in- 
controvertible. 6 - ) Theorists have also predicted that, un- 
like weak topological insulators, an STI exhibits a sin- 
gle surface Dirac cone irrespective of the orientation of 
the surface. However, the ARPES measurements can be 
carried out only on the perfectly cleaved surface of lay- 
ered STI samples, i.e., ARPES is restricted to surfaces 
of some particular easy-to-cleave orientations. It is in- 
deed a much less trivial issue whether the gapless Dirac 
state observed on such one "good" facet of a crystal is 
smoothly extended to adjacent ones, eventually covering 
the entire surface of the sample. Experimentally, such a 
behavior of the surface Dirac state may be accessible by 
STM measurements. In ref. 7, an example of such an ex- 
periment performed on atomic-scale terraces of an STI 
has been reported. 

The protected gapless surface Dirac state is known to 
be robust against perturbations that do not break time- 
reversal symmetry. An electron in a protected surface 
state exhibits a notable feature, that is, its spin is locked 
to a particular direction determined by its momentum 
(spin-to- momentum locking). This implies that the spin 
direction of an electronic state with wave number k is or- 
thogonal to that with wave number — fc, resulting in the 
complete suppression of 180° backward scattering. Here 
we concern ourselves with the transmission through an 
interface of two Dirac electron systems, 8 ' in which the 
absence of 180° backward scattering naturally plays a 
central role. Generally, gapless Dirac electrons, such as 
those found in graphene, are not necessarily immune to 



backward scattering at the interface 9, 10 - ) if the incom- 
ing electron is away from the normal incidence. Only an 
electron normally incident to the interface is forced to be 
completely transmitted by the absence of 180° backward 
scattering. The greater the deviation from the normal 
incidence, the larger the reflection becomes. 




Fig. 1. Concave 90° step edge consisting of two flat surfaces of 
a 3D STI: one surface is on the yz plane, and the other on the zx 
plane, and they meet at the 2-axis. 

Let us focus on the behavior of the protected sur- 
face state on the L-shaped wall of an STI (see Fig. 1). 
This structure can be regarded as a concave 90° step 
edge consisting of horizontal and vertical facets. Note 
that both the horizontal terrace and the vertical wall 
exhibit a protected gapless Dirac cone, since they both 
separate an STI and a vacuum, two topologically dis- 
tinct worlds. Although an STI exhibits a protected gap- 
less Dirac cone on surfaces of an arbitrary orientation, it 
can have, and in practice, it always exhibits anisotropy 
in its model parameters. As a natural consequence of 
this, Dirac cones on facets of different orientations (here, 
those on the vertical and horizontal surfaces) generally 
have different apertures exhibiting different group veloc- 
ities. With this taken into account, are these two Dirac 
cones still smoothly connected without being affected by 
scattering at the interface, even away from the normal 
incidence? Naturally, this is an issue closely related to 
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the corresponding problem in the junction of two pure 
2D Dirac systems such as in graphene. In the case of the 
bulk-surface system we consider, is there any difference 
in the immunity of the Dirac electrons from scattering 
processes? These are the issues we would like to address 
in this letter. We set h = 1 below. 

To state the problem unambiguously, let us first con- 
sider the surface electronic states in the asymptotic re- 
gions sufficiently away from the step edge. We focus on 
a single concave 90° step edge consisting of two semi- 
infinite surfaces, one on the yz plane and the other on the 
zx plane (see Fig. 1). The system is assumed to be trans- 
lationally invariant in the z-direction. By treating each 
asymptotic region individually, the effective Hamiltoni- 
ans on the yz and zx planes can be respectively deduced 
as 
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We take into account the anisotropy of the (bulk) sys- 
tem, which results in the velocity mismatch, i.e., if 
A x 7^ Ay. We focus on cigcnstates with energy E = 



[{Akf + (Bk,)} 1 / 2 , where Ak ; 
responding eigenfunctions are 



Ayky 



A x k x . The cor- 
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where e ±lx = {Ak ± iBk z )/E and the prefactor is intro- 
duced to normalize the probability current in the direc- 
tion perpendicular to the z-axis. Here, the sign ± indi- 
cates the direction of propagation (see Fig. 1). Let us 
consider the scattering problem in which is incident 
from y = oo. If k z — 0, the reflection into a^ z is forbid- 
den owing to the absence of 180° backward scattering, 
and hence the incident wave is completely transmitted 
to a?j?. Here comes the question to be answered: n ) does 
reflection occur when k z ^ 0? 

Note that solving a scattering problem at a junction 
of two inequivalent (velocity-mismatched) Dirac systems 
is a formidable task. Since one of our primary interests 
is to examine the continuity of the current density at 
the corner, we need to deal with a subtle issue related 
to the discontinuity of the wave functions, which is fre- 
quently encountered in the study of Dirac systems. 8 ~ 10 ) 
To overcome such a difficulty, we introduce the following 
"hyperbolic model" (the details of the model specified 
later), in which the surface of our 3D STI system with 
a 90° concave step edge is regarded as a limiting form 
of a hyperbola (Fig. 2). The two asymptotes of the hy- 
perbolic surface represent the horizontal terrace and the 
vertical wall meeting at the concave right angle of the STI 
crystal. Using this hyperbolic model and a set of contin- 
uous curvilinear coordinates in terms of this hyperbola, 
we derive the surface effective Hamiltonian valid over 
the entire surface, i.e., this provides us with a scheme 
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Fig. 2. Cross section (on the xy plane) of our "hyperbolic 
model" , representing a curved surface of a 3D STI. The cross sec- 
tion of the hyperbolic surface is given by xy = R 2 , reproducing a 
90° step edge in the limit of sufficiently small R. 



for treating the two asymptotic surfaces in the frame- 
work of a single Hamiltonian. 13 ) Note that in the con- 
ventional approach, surface electronic states in the two 
asymptotic regimes, such as those given by eqs. (3) and 
(4) , are derived separately as a solution of two individual 
uncorrelated problems. Here, we treat them on the same 
footing to discuss (eventually) how they are connected 
at the step edge. The treatment of the opposite 90° con- 
vex step edge introduces further technical complications 
in the proposed approach. This last issue will be briefly 
discussed toward the end of the paper. 

Let us start with the following bulk effective Hamilto- 
nian for 3D anisotropic STIs in the continuum limit: 14, 15 ) 
#buik = m(p)T z + (A x p x a x + AyPyOy + Bp z a z )r x , where 
p, = -id, (i = x,y,z) and m(p) = m + m 2 (pl + p 2 y + p 2 z ) 
is the mass term. Without loss of generality, we assume 
that mo > and iri2 < 0. The two types of Pauli ma- 
trices er = (cr x , <Ty, <j z ) and r = {t Xi t Vi t z ) respectively 
represent the real and orbital spin degrees of freedom. If 
the ordinary matrix representation of <x is used, i?buik is 
expressed as 



#bulk 



m{p)T z + Bp z r x 
V+t x 



V-T X 

m(p)T z - Bp z r x 



, (5) 



where T>± = A x p x ± iA y p y . It should be noted that H^g 
and "H z | can be derived from -ffbuik- 14, 15 - ) 

With the basic assumptions stated, let us introduce 
an ingenious trick. As mentioned earlier, the main point 
of our protocol is to consider, instead of directly treat- 
ing the 90° step edge, the hyperbolic system depicted in 
Fig. 2. We assume that the intersection of the curved 
surface with the xy plane is a rectangular hyperbola, 
xy = R 2 , and that the 3D STI is translationally invariant 
in the z-direction. If R is sufficiently small, this surface 
can be regarded as a concave 90° step edge. It is conve- 
nient to introduce curvilinear coordinates, as shown in 
Fig. 2. Let us draw a straight line that perpendicularly 
crosses the hyperbola at the crossing point (xo,yo)- We 
define <j> as the angle between the x-axis and this line, and 
r as the distance from the crossing point. The Cartesian 
coordinates (x, y) are expressed as 



x = —rcoscf) + Xo(<p), 



(6) 
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y = -r sin <j) + y ((f)), 



(7) malization constant, and 



where x = i?tan5 and y = Root? <\>. The infinitesi- 
mal cross-sectional area is given by dS = (r + f ((/>)) drd(f>, 
where 

R 



2sm 2 0cos2 



The derivatives are expressed as 

d x — — cos <j)d r + (r + sin 



■ sin (j)d r — (r + /) 1 cos 



(8) 

(9) 
(10) 



With these expressions, we can rewrite the mass term as 
m(p) = m± + mil with 

m ± = m -m 2 [d 2 . + (r + fy 1 d r ] , (11) 

m 2 [(r + f)- 2 dl-(r + f)- 3 (d rt> f)d lt> ]. (12) 



mil 



Similarly, V± is rewritten as T>± 
where 

n ± =i(r + f)- 1 (-A x sin<j)±iA y cos(f>)d 4n (13) 

e^ = (Ascos^ii^sin^)/^ (14) 

with A^ = (Al cos 2 <p + A 2 y sin 2 <j)) 1/2 . 

To derive the surface effective Hamiltonian in the spirit 
of k ■ ^-approximation, we divide i?b u ik into two parts as 
#buik = H±_ + H 11 , where 



Hs = 



m±T z iA<pe ^d r T x 



m»T z + Bp z r x 



n + r x 



m»T z - Bp z r x 



(15) 



(16) 



We first solve the radial eigenvalue equation 14-17 ) 
H±\tp) = Ej_\ip) with the boundary condition \ip(r = 
0)) = 0, that is, all four components of the wave func- 
tion |?/>) vanish on the surface. As the simplest approxi- 
mation, we replace (r + f) d r in m± with ((r + f)~ 1 )d r , 
where the definition of the average ((r + is given 

below [see eq. (26)]. Then, we can show that the eigen- 
value equation has surface solutions of the damped form, 
|V) = c~ Kr |w), where k _1 measures the penetration of 
the surface wave function into the bulk. By superposing 
two damped solutions, 18 ^ we construct the solution of 
the radial eigenvalue equation localized near the surface 
as l^) = c~ K - r \u-) — c _K+r |w + ). The boundary condi- 
tion |\I>(r = 0)) = holds only when = \u + ) for 
k + 7^ k_. As shown in ref. 17, this results in the zero- 
energy condition E± = in our model. We thus find that 
k± is given by 



k±(4) = (A4, ± yjAl + 4m m 2 ) /(-2m 2 ) (17) 

with Atj, = A^ — m 2 ((r + We also find that two 

basis eigenstates, |+)) and |— )), for H± with zero eigen- 
value are given by 

1 



l±» 



-_p{r, 4>)\n±), 



(18) 



where p(r, </>) = e 



- c~ K+r , is a ^-dependent nor- 



\n±) = - 



e -i*/2 
Te i0/2 



1 

Ti 
1 



(19) 



The real-spin sector \S±) = '(e-^ 2 ,^'^ 2 )/^ points 
in the ±s x direction in the regime of yo — > 00 (<f> — > 0), 
while it points in the ±s y direction in the regime of xq — > 
00 ((f) — > it/ 2). It should be emphasized that when R is 
very small, \S±) rapidly changes its direction from ±s x 
to ±s y across the point (xq, yo) = (R,R). 

Within the k -p-approximation, any surface state \ip}} 
can be represented as a linear combination of |+)) and 
|— )) with the amplitude respectively specified by a+ 
and a_, i.e., = a + |+)) + a_|— )). The effective 

surface Hamiltonian "H c ff for the two-component spinor 
a. = *(a + ,a_) is defined by 



H eS = 



((+\H\\\+)) 



«-|^|||-» 



Here, each matrix element is expressed by 

p(r, 4>) 



<W|I*'» 



dr(r + /) 



na\Hu\na') 



(20) 



p(r, <p) 



c<p 
(21) 

where the factor r + f reflects the fact that dS = 
(r + f)drd(j) and = J °° dr(r + f)p(r,(f>) 2 . We easily 
find that ((±|ff|||±)) = 0. In evaluating the off-diagonal 
elements, we should note that d^, in H\\ acts on not only 
|ner') but also p(r, 4>)/^c^. After tedious but straightfor- 
ward calculations, we find 



H eS = 



where 



«±I#|||T» 



with 



((+\H\\ 

«-|#nl+» 



(22) 



A 



(r) + f 



(r)+f 



±iB Pz 



A{<j>)=[A 4> -m 2 {{r + f)- 1 )] 
In these equations, we have used the notation 
Jo°° drrp(r, (f>) 2 



(r) 



Jo°° dr P( r > ' 



((r + fr 1 ) = 



J Q °°dr(r + f) V(r,i 



(23) 



(24) 



(25) 



(26) 



The second term on the right-hand side of cq. (23) is 
essential in ensuring the hermiticity of % e ff- 

It is convenient to introduce the one-dimensional co- 
ordinate 



r<t> 

= / ( 

J-k/4 



<¥Kr>(0') + W)]. 



(27) 



located slightly beneath the (geometrical) surface (see 
Fig. 2). Note that the limit of 4> — > (n/2) corresponds 
to I — > —00 (+00). We rewrite the effective Hamiltonian 
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with this coordinate. Since dl = ((r) + f)d(f>, the wave 
function a (I) in the new coordinate is related to the old 
one by a(l) = a{4>)/{{r) + f) 1 / 2 . From eq. (22), we find 
that the effective Hamiltonian for a (I) is expressed as 



H 



eff 



A P i 





\ (vi~A) 



A P i + 



3 M 



iB Pz 



iBp z 



(28) 

with pi = —idi. This form of an effective Dirac Hamilto- 
nian has been suggested in ref. 8. The probability current 
operator in the I direction is given by 



Ji = 



A 
A 



(29) 



Let us consider the behavior of A(l) given in cq. (24), 
which should be regarded as the effective velocity along 
the one-dimensional coordinate. Since / — > oo in the 
limit of I — > ±oo [sec cq. (8) and the note given below 
cq. (27)], we observe with the aid of d^cj) = A x A y jA\ 
that A(— oo) = A y and A(po) = A x . Contrastingly, the 
velocity can become very large in the close vicinity of 
the step edge if R is vanishingly small. Note that m 2 < 
is assumed from the outset, and that / « near I = 
(i.e., 4> = 7r/4) if R is very small. The correction term 
— 1712 ((r + m A becomes very large, leading to a 

large increase in velocity. 

Now, we can answer the main question raised earlier. 
The solution of H e sot = Ea with E = [(Ak) 2 + (Bk z )} 1 / 2 
is obtained as 



1 



a± = 



2A(l) 



1 

±e Tix 



exp ±i 



i 

dl'4 
o A(l') 




We see that eq. (30) respectively reproduces eqs. (3) and 
(4) in the regimes of / <C —R and I 3> R. We also see from 
cq. (29) that the prefactor A(l) -1 / 2 guarantees current 
conservation over the entire system. That is, this ct± con- 
tinuously connects the two asymptotic eigenfunctions. 19 ) 
In addition, eq. (30) is justified even for vanishingly small 
(but finite) R. Therefore, we conclude that no reflection 
takes place at a concave 90° step edge, although the ve- 
locity becomes very large in its vicinity and hence the 
amplitude of the wave function is reduced. 

Let us now consider an inverted situation in which 
the convex side of the hyperbola (see Fig. 2) is filled 
with a 3D STL In this situation, our curvilinear coor- 
dinates (r < and < <j> < 7r/2) can be applied only 
to a limited region near the hyperbola of width much 
smaller than R. Indeed, d x and d y in eqs. (9) and (10) 
become ill-defined owing to the presence of (— |r| + /) 
if M ^ R- Therefore, we restrict our analysis on sur- 
face states to the case where R is sufficiently larger than 
the penetration depth (|r|). Analysis similar to that re- 
ported above reveals that surface states in the convex 
case are described by an effective Hamiltonian essentially 
equivalent to eq. (28) but the velocity, now given by 
A(4>) = [A$ + m 2 ((-\r\ + d^, is reduced [note 



the sign change in front of the correction term compared 
with that in eq. (24)]. This implies that the behavior of 
surface states in the convex case is essentially equivalent 
to that in the concave case except that the velocity is 
renormalizcd in the opposite way. 

We finally comment on the validity of our analysis. 
Since -ffbuik is valid in the long-wavelength regime, one 
may question whether it can be applied to STIs with a 
sharply edged surface. In eq. (30), the shortest length 
scale of the derived eigenfunction is on the order of (r) 
even though R becomes vanishingly small. Thus, we ex- 
pect that our approach will be justified as long as (r) is 
much longer than the lattice constant. This does not nec- 
essarily mean that the employed approach allows quan- 
titative predictions of the electronic properties of, for ex- 
ample, atomic-scale terraces of an STI with atomic-scale 
precision. Yet, our conclusion (no reflection at a single 
step edge) itself is consistent with the experimental result 
indicating that topological surface states are transmitted 
through atomic-scale steps with high probabilities. 7 ) 
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